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We present conservative 3+1 general relativistic variable Eddington tensor radiation transport 
equations, including greater elaboration of the momentum space divergence (that is, the energy 
derivative term) than in previous work. These equations are intended for use in simulations involving 
numerical relativity, particularly in the absence of spherical symmetry. The independent variables 
are the lab frame coordinate basis spacetime position coordinates and the particle energy measured 
in the comoving frame. With an eye towards astrophysical applications — such as core-collapse 
supernovae and compact object mergers — in which the fluid includes nuclei and/or nuclear matter 
at finite temperature, and in which the transported particles are neutrinos, we pay special attention 
to the consistency of four-momentum and lepton number exchange between neutrinos and the fluid, 
showing the term-by-term cancellations that must occur for this consistency to be achieved. 
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I. INTRODUCTION 



Neutrino transport is a necessary ingredient of core- 
collapse supernova simulations [l|-|3|. Determining the 
fate of the stellar material — for instance, does an explo- 
sion happen, and if so, how? — requires calculation of the 
four-momentum and lepton number exchange between 
the fluid (which includes nuclei and / or nuclear matter at 
finite temperature) and the neutrinos that stream from 
and through it. For the purpose of studying the explo- 
sion mechanism, we take the traditional approach and 
consider only massless neutrinos described by classical 
distribution functions (phase space densities) /(i,x, p). 
(Calculation of the emerging neutrino signals — of intrin- 
sic interest as an observational probe of the core-collapse 
supernova environment, and of the properties of the neu- 
trinos themselves — definitely requires treatment of the 
quantum effects induced by neutrino mass and flavor 
mixing Recent explorations suggest that flavor 

mixing does not impact the explosion mechansim @-[l3| , 
but consensus on the impacts of flavor mixing in super- 
novae has a fickle history, and future more definitive sim- 
ulations that include neutrino transport with quantum 
kinetics could surprise us with flavor mixing effects on 
the explosion mechanism as well.) 

The formidable challenge of solving for the neutrino 
distribution functions /(i,x, p) can be ameliorated by 
solving only for their angular moments, with some form 
of closure serving to truncate the scheme at low or- 
der. Common examples include flux-limited diffusion 
14], 15]. 16], truncated at the zeroth moment J(t,x,c). 
with prescriptions for the first and second moments 
% l (t,yL, e) and JC^jt, x, e); and a variable Eddington ten- 
sor approach [17H23| . truncated at the first moment 
ri l (t,x,e), with prescriptions for the second and third 
moments K}\t, x, e) and C Vlk {t, x, e) (see Sec- ITXC]) . (Note 



that Refs. [19|, |21j, |22| integrate out the neutrino energy 
as well as angle dependence, which entails several ma- 
jor simplifications: in addition to the further reduction 
of dimensionality, it eliminates neutrino-energy-changing 
interactions, velocity-dependent terms, and the third mo- 
ment C 1 ^ . On the importance of the neutrino-energy- 
changing interactions and velocity-dependent terms, see 
Ref. gj) 

Notably for the focus of this paper, conservative 3+1 
general relativistic variable Eddington tensor radiation 
moment equations are presented by Shibata et al. [20j . 
Their variables are functions of lab frame coordinate ba- 
sis spacetime position coordinates a; M — that is, the space- 
time coordinates that appear in the 3+1 metric — but the 
momentum space dependence is on the neutrino energy 
e measured in the comoving frame. The long-recognized 
freedom to choose different coordinate systems f or sp ace- 
time and momentum space (e.g. Refs. [25M3l| ) al- 
lows particle/fluid interactions to be evaluated in the co- 
moving frame in the context of Eulerian grid-based ap- 
proaches to multidimensional spatial dependence. 

Apart from the fact that we leave aside for now dis- 
cussion of specific closure prescriptions and contributions 
to the collision integral (which are stories unto them- 
selves) , one difference between our presentation from that 
of Shibata et al. is our starting point. Like En- 

deve et al. [23l | . we begin from conservative reformu- 
latio ns |32| of th e g eneral relativistic Boltzmann equa- 
tion [25l |29|. [30l [33 . |34| rather than the moments for- 
malism of Thorne 35|- By angular integration of the 
four-momentum conservative reformulation of the Boltz- 
mann equation, we obtain in Sec.|TT]a general relativistic 
variable Eddington tensor formalism in which the rela- 
tionship between the lab frame (denoted by unadorned 
indices) and the comoving frame (denoted by hatted in- 
dices) is expressed in terms of coordinate transformations 
L^k and comoving frame connection coefficients T^pf, 
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(see also Ref. H|). 

In specializing to the 3+1 formulation of general rela- 
tivity in Sec. we extend the treatment of Shibata et 
al. [20l | with a full elaboration of the momentum space 
divergence (i.e. the energy derivative term, in this angle- 
integrated moments case). Important aspects of our ap- 
proach include (a) consistent use of what we call 'Eu- 
lerian decompositions' and 'Eulerian projections,' which 
are natural to the 3+1 approach; and relatedly, (b) a 
shift from conceptualizing the relationship between the 
lab and comoving frames from coordinate transforma- 
tions to the (covariant) relative three-velocity 
connecting the four- velocities n M and of Eulerian and 
Lagrangian observers. Our approach in Sec. IIIII is more 
geometric than that in Sec. [TT] (in conception if not nota- 
tion); indeed it allows us to obtain explicit results while 
almost completely avoiding encounters with connection 
coefficients. 

We also add to the treatment by Shibata et al. [2(| 
by showing, in both Sec. [TT] and Sec. IIIII how the four- 
momentum exchange with the fluid expressed by a con- 
servative variable Eddington tensor formalism is consis- 
tent with a conservative treatment of lepton number ex- 
change (see also Ref. [23| ). A conservative treatment of 
four-momentum exchange with the fluid, properly dis- 
cretized for consistency with conservative number ex- 
change, is expected to facilitate simultaneous energy and 
lepton conservation in numerical simulations — an impor- 
tant check on the physical reliability of simulation out- 
comes [3(| . In this respect, there may be room for 
improvement over cases in which this consistency has 
not been considered (e.g. Refs. [U HH), or in which 
consistency between conservative number exchange and 
non-conservative four-momentum exchange has been ad- 
dressed (e.g. Ref. [18]). The nature of the consistency 
of our conservative four-momentum transport equations 
(modulo gravitational sources) with a conservative num- 
ber transport equation is made particularly explicit in 
Sec. IIIII in which we elucidate the term- by-term cancella- 
tions that must occur for this consistency to be achieved. 



II. GENERAL RELATIVATISTIC VARIABLE 
EDDINGTON TENSOR FORMALISM 

After exhibiting the Boltzmann equation and its con- 
servative reformulations in terms of number and four- 
momentum exchange with the fluid, we obtain from 
the latter a variable Eddington tensor formalism in 
which the relationship between the lab frame (denoted 
by unadorned indices) and the comoving frame (de- 
noted by hatted indices) is expressed in terms of coor- 
dinate transformations L^^ and comoving frame con- 
nection coefficients T^^p, and show how it relates to 
the angle-integrated number-conservative reformulation. 
This ground is also covered in Ref. [23| , but is considered 
here to establish notation and provide a self-contained 
presentation. 



A. The Boltzmann equation 

Classical neutrino distribution functions f(t, x, p) are 
governed by the Boltzmann equation [25|, [29|, [30, [34| . 
In its geometric form, it states that the change in / along 
a phase space trajectory with affine parameter A is equal 
to the phase space density C[f) of point-like collisions 
that add or remove particles from the trajectory: 



dX 



C[f]- 



(1) 



The phase space measure is defined in such a way that / 
and C[f] are both invariant scalars. For practical com- 
putations it is necessary to introduce phase space coordi- 
nates: spacetime coordinates x M , and momentum space 
coordinates p 1 (the timelike momentum component p° is 
fixed in terms of the spacelike components p l by the mass 
shell constraint). In terms of these coordinates, Eq. (fTJ) 
becomes 



efa+ df dp 1 df 
~dXdx^ + 'dXdp 1 



= c[f}. 



The geodesic equations describing the trajectory are 

~dX 

dp» 

~dX 



if— = _r M v v v p 



so that Eq. ([2]) becomes 



V 



df_ 

<9a+ 



(2) 

(3) 
(4) 

(5) 



now an integro-partial differential equation (the integrals 
appearing on the right-hand side). 

There is freedom in choosing the spacetime and mo- 
mentum space coordinates. Taking the unadorned in- 
dices to denote what we shall call a lab frame coordi- 
nate basis (also called a 'natural' or 'holonomic' basis), 
the connection coefficients T^ p are given in terms of the 
spacetime metric g pv as 



r 



vp 9 .y 



f dg<jv dg ap _ dg Vl 
\ dx p dx v dx° 



(6) 



However, it is most convenient to express the particle in- 
teractions entering C[f] in terms of momentum compo- 
nents p % reckoned with respect to an orthonormal refer- 
ence frame comoving with the fluid (a 'comoving frame'). 
We define a composite transformation 



L I* — 



A A 1 . 



(7) 



consisting of a Lorentz boost A^ from an orthonormal 
comoving frame (denoted by indices with a hat) to an 
orthonormal lab frame (denoted by indices with a bar), 
followed by a transformation to the lab frame coordinate 
basis with a local tetrad e^p. This tetrad is independent 
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of the fluid velocity; it locally transforms the metric into 
the Lorentz form (r^p) = diag[— 1, 1,1,1]: 



fl V 



(8) 



Of course, the boost preserves the Lorentz metric; 
this implies that the composite transformation is 
itself also a tetrad. The inverse of Eq. is 



(9) 



expressed in terms of the inverse tetrad e' 3 M and inverse 
boost A^. In terms of lab frame coordinate basis space- 
time components and comoving frame momentum com- 
ponents, the Boltzmann equation reads 

L,1 A/^r-r^//^ = c[/] ! (io) 

where the connection coefficients in the comoving frame 
are 



(11) 



Finally, assuming particles of zero mass, it is convenient 
to express the comoving frame null momentum compo- 
nents in terms of energy, polar angle, and azimuthal 
angle, that is, in terms of momentum space spheri- 
cal polar coordinates (denoted by indices with a tilde) 
(p x ) = M, iff: 



(/)=e(MVV § ) : 



= e(l, cos sin $ cosy, sin § sin (p) , (12) 

which also defines the unit normal three- vector t tangent 
to the comoving- frame three- momentum p % . In terms of 
these momentum space coordinates the Boltzmann equa- 
tion now reads 

- df - - -<V~ df 

where 

—j — — — sind cos $ cos tp cos ■& sin ip (14) 
P e \ — sin <p/ sin?? cos ip/ sin i) I 

is the Jacobian relating momentum space spherical and 
Cartesian coordinates. 



B. Conservative reformulations of the Boltzmann 
equation 

Conservative reformulations of the Boltzmann equa- 
tion are available 32j that render plain its connection to 
number and four-momentum conservation (or balance, 
given the presence of source terms), and therefore may 



be helpful in attempts to maintain fidelity to global con- 
servation laws in numerical simulations. 

The number-conservative reformulation of Eq. (fi"3"j) is 



>N 



M N = C[f], 



(15) 



with spacetime divergence Sn and momentum space di 
vergence Mjv given by 

1 d 



Sn — 
M N = 



l—g dx^ 1 
1 d 



. , . . -esintfr^^pVV 
esinz? dpi \ Ufi dp 1 



(16) 



(17) 



When Eq. (|15p is integrated over the invariant momen- 
tum space volume element (e.g. Ref. [251 ] ) 

, D i — dip* d 2 p j d 3 p k 

dP = V-9£tjk 7 : (18) 

K-po) 

= esin-dded-ddip, (19) 

the momentum space divergence term manifestly disap- 
pears, leaving the number balance equation 



1 
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-g dx^ 



(V=gN») = / c[f 



dP 



where 



N" = / p»f 



dP 
(2nf 



(20) 



(21) 



is the number flux vector (e.g. Ref. [25|), expressed here 
in terms of the lab frame coordinate basis (note p^ — 
L^'fip' 1 )- We use units in which h = c = 1; relative to 
works in which instead h — c — 1 , this leads to the factors 
of (27r) 3 in the preceding two equations. 

Similarly, the four-momentum-conservative reformula- 
tion of Eq. ([T3j is 

(S T ) p + (M T ) p = L" pI fC[f], (22) 

with spacetime divergence St and momentum space di- 
vergence Mt given by 

(S T y = ^=^77 (V^gLPpL^ptp* /) 



(M T y 



/—g dx^ 

r^L^L^ptptf 

1 d 



(23) 



e sin d dp 1 



dp> 



(24) 



When integrated over dP, Eq. ([22]) yields the four- 
momentum balance equation 

1 f) f dP 

(J—gT^) + Kp t»» = j P " c[f] — ^ , 



-g dx^ 



where 



dP 
(2nf 



(25) 



(26) 



is the stress-energy tensor (e.g. Ref. |25|). expressed here 
in terms of the lab frame coordinate basis. 
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C. Variable Eddington tensor formalism 

Solving for / in its full dimensionality being computa- 
tionally overwhelming, the dimensionality of the problem 
can be reduced by considering only its lowest angular mo- 
ments. A truncation of the hierarchy of moments must 
be performed by means of closure relations (see e.g. the 
references in Sec. HJ. Just as it is most convenient to de- 
scribe neutrino interactions with the fluid in terms of mo- 
mentum components reckoned in the comoving frame, so 
also it seems sensible to define angular moments and pre- 
scribe closure relations in the comoving frame. Wc define 
the lowest angular moments of / (x M , e, n) as follows: 



J(x^e) = e, 


f f(x^e,a) dn, 


(27) 


W{x^e)=e 


f tf{x^e,si) dn, 


(28) 


/C ?J "(y\e) = e j 


f f{x»,e,fL) dfl, 


(29) 


C^(x^e) = e J 


f tik k f(x»,e,n) dn, 


(30) 



where l % is the unit three-vector tangent to the three- 
momentum in the comoving frame, defined in connection 
with Eq. (|12l) . The integration over dn = sin t9 dd d<p is 
performed over the unit sphere. Note that the energy 
dependence is retained; these monochromatic moments 
are functions of lab frame coordinate basis spacetime po- 
sition components x^ and the comoving frame energy e. 
(This is the first instance of a convention we employ, of 
denoting monochromatic or energy-dependent quantities 
with script symbols.) The flux-limited diffusion approx- 
imation entails truncation of the hierarchy at the zeroth 
moment J , with prescriptions for the first and second 
moments H z and K? J in terms of J . In the variable Ed- 
dington tensor approach the hierarchy is truncated at TV , 
with the next higher moments rewritten as 



that is, 



ffk 



}C ] = k t: >J, 

, _ J m$ f dn 
= j f dn ' 

J tiH k j dn 



(31) 
(32) 

(33) 
(34) 



jfdn ■ 

A number of different approaches to computing the Ed- 
dington tensors k v and l lJ ~ k mig ht be taken; see for in- 
stance the references in Ref. [23j j . 

We choose the variable Eddington tensor approach, 
and note that we require four equations for the four un- 
knowns J, W (in addition to whatever scheme is used to 
compute the Eddington tensors). Inspection of Eqs. (|22)l - 
(|26p indicates that the four-momentum conservative for- 
mulation of the Boltzmann equation may serve as a suit- 
able basis for the four equations we require. They are 



suggestive of conservative evolution of zeroth and first 
moments in the lab frame coordinate basis, which may 
prove helpful in maintaining numerical four-momentum 
conservation; yet the arguments of the spacetime and mo- 
mentum space divergences can nevertheless be expressed 
in terms of the comoving-frame moments J , "H 1 we must 
take as our primitive unknowns. 

We prepare to implement this strategy with some ad- 
ditional definitions. We define a monochromatic stress 
energy, whose components are functions of lab frame co- 
ordinate basis spacetime position components x* 1 and the 
comoving frame energy e: 



pV7 fy\ e, n) dn. (35) 



Its components are related to the comoving frame mo- 
ments by 




J W 



J W 
W 0J 



Similarly we define 



U 



(x'V) - 



/pV/^.e.fl) dn, 



(36) 



(37) 



whose components are given by 

U 6 ^ = = W ob =eT^, (38) 
u m = e£ «fe = e ffkj ( 39 ) 

The bottom line is that all the components of both 
and W >,vp are just our primitive unknowns J and % l , 
modulo factors (taken to be known) of e, k^ , and l %3 . 

We obtain our equations for J and H l by integrating 
Eq. (|2"2"|) over dn and dividing by e: 



(S T ) P + (M T )" = L p f 



c[f] dn, 



(40) 



where the angle-integrated spacetime divergence (Sy) p 
and momentum space divergence (Mr) p (now denoted 
in a sans-serif font to distinguish them from the unintc- 
grated (St) p and (Mt) p ) are given by 



(St)P = 7^^^ l ^ t ' a ) 
(M T ) P = 1| Lev, J ri^pWfdn 



To further simplify (Mj) p , note that 

dp 1 _ Pj_ 
dp 1 e 



(41) 



(42) 



(43) 
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by virtue of Eq. (1121) and the first row of Eq. (1141) : and 
that 11 III 



p/i PiP P 



(44) 



This equation, which follows from = d(pi- l pf^J /dX = 
2pf l dpi 1 1 dX and Eq. ((3]), is important because it makes 
only the third (rather than fourth) moment appear [23j . 
With these relations we obtain 



(M r ) 



" - — — (-e 2 L p - r? - 
e 2 de { e L P l ^ U 



(45) 



for the angle-integrated momentum space divergence. 

Equations (|40|) . (|4Tj) . and (|45]) are the relations we 
seek. They provide four equations (p = 0,1,2,3) for 
our four primitive unknowns J and TV, the comoving 
frame angular moments; this is because the components 
of T^ v and W" /p are in fact J and W l , modulo factors 
(taken to be known) of the comoving energy e and Ed- 

dington tensors k 1 ' 1 and (see Eqs. ([551). . and (j3^ |) ) . 
Moreover, Eqs. (|40l) . (l4Tj) , and (1451) are conservative in 
(that is, expressed in terms of divergences with respect 
to) our chosen phase space coordinates — the lab frame 
coordinate basis spacetime position coordinates and 
the comoving frame neutrino energy e. 



D. Four-momentum and lepton number exchange 

The right-hand side of Eq. (|4T))) is intimately related 
to source terms for the fluid energy and momentum 
equations. From the right-hand side of Eq. (|40[) . the 
monochromatic four-momentum source (f a for neutrino 
species a is 



Qa 



pi' 



c a [f] dn. 



(46) 



Integrating e 2 /(27r) 3 times Eqs. gO), (J4l]), and (J45]) over 
neutrino energy e for a particular neutrino species a, we 
have 



V7 rpv\l 



e 2 de 
(2tt)3 



(47) 



where V M denotes the covariant derivative. The diver- 
gence of the total stress-energy — which includes the fluid 
and all species of neutrinos — must vanish: 



By virtue of Eq. (|4T|) we have 



(48) 



(49) 



which expresses the four-momentum exchange between 
the fluid and the neutrinos. 



To address the exchange of electron lepton number we 
must consider the number-conservative neutrino equa- 
tion. We begin by defining a monochromatic number 
flux 



A^(:r/\e) 



//(xV,^ dn. 



(50) 



Inspection of Eqs. (fT2"j). ([331) and (|3"o]) shows that it is 
related to the monochromatic stress energy and comoving 
moments by 

(/v*) = I (r 5 ^) = - (t^) = - (j, wf . (5i) 



Following steps similar to those used to obtain Eqs. (|40[) . 
(l4"TT) . and (|4"5j) . we integrate the number conservative 
Eq. (p~5|) over cf£l and divide by e to obtain 



- / c[f] dn, 



(52) 



where the spacetime divergence S n and momentum space 
divergence Mn are given by 



1 



d 



5jV 



I AT 



-<7 
l__d / 



2. pU n^y^ 



(53) 



(54) 



We define a monochromatic lepton number source r a for 
neutrino species a from the right-hand side of Eq. (|52[) : 



c„[/] dn. 



(55) 



Integrating e 2 /(27r) 3 times Eqs. ([52"j) - (f54|) over neutrino 
energy e for neutrino species a, we have 



(2nY 



= Ra- 



(56) 



The divergence of the total electron lepton number van- 
ishes: 



V M (N? + iV£ - JV£) = 0, 



(57) 



where iV** is the net electron number flux vector of the 
fluid. By virtue of Eq. (f5"6"j) we have 



(58) 



which expresses electron lepton number exchange be- 
tween the fluid and the neutrinos. 

Moreover, the lepton number source on the right- 
hand side of Eq. (|S"5|) is intimately related to the four- 
momentum sources. Comparing Eqs. (|55[) and (|46l) . we 
see that 



e 



(59) 



(note that iP^IM-n = <5%). 
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The numerical consistency of a scheme in which the 
lepton number exchange is expressed in terms of the four- 
momentum exchange via Eq. ([59)) — which consistency 
presumably has consequences for simultaneous conserva- 
tion of four-momentum and lepton number [36j ] — depends 
on the discretizations chosen for Eqs. (j40|) . (|4T]) . and (|45|) . 
(This is true regardless of whether the fluid/neutrino 
couplings are handled in an operator split fashion, e.g. 
Refs. [13, E, H3|; or simultaneously with the solution 
of the transport equations in a single implicit solve, e.g. 
Refs. [3^,[3g].) To see this, note that the relation em- 
ployed in Eq. (151?)) between the right-hand sides of the 
number and four-momentum equations should apply to 
the left-hand sides as well. In particular, the identity (see 
also Ref. (H) 

S N + M N = h (> p [(S T ) p + (M T ) p } (60) 

should hold for the discretized equations. Examining the 
first term on the right-hand side, we find using Eqs. (pTTj) 
and EB that 



-L 6 P (s T r 



L n L u L»,,rP 




~ LP ° L "^dx» 



7* A . (61) 



For the second term on the right-hand side of Eq. (f60|) . 
we find using Eqs. (g5J and d3HJ) that 



L° p (M T r 



M 



(62) 



Thus we see that in the sum of the above two equations, 
the 'extra' terms on the right do indeed cancel by virtue 
of Eq. (HI]). (Note that L p <9 M L 6 p = -L 6 p d p L p , thanks 

to = d p (<5%) = d p (iPpLPj^ .) Ideally, the analytic 

steps confirming Eq. (l60l) can be followed in the discrete 
limit in order to find a discretization of Eqs. (|40|) . (|4T]) . 
and (j45|) that is consistent with respect to both four- 
momentum and lepton number exchange. 



III. SPECIALIZATION TO THE 3+1 METRIC 

After briefly reviewing the 3+1 formulation of gen- 
eral relativity, we show how thinking in terms of three 
vectors — the four velocity n u of Eulerian observers, the 
four-velocity u M of Lagrangian observers, and the (co- 
variant) relative three-velocity that connects them — 
facilitates a full elaboration of the 3+1 variable Edding- 
ton tensor moment equations, including the detailed re- 
lationship between these and the number exchange equa- 
tion. 



A. Description and evolution of the geometry 

Numerical relativity often is built upon the 3+1 formu- 
lation of general relativity. In this approach one considers 
a foliation of spacetime into spacelike slices, i.e. three- 
dimensional hypersurfaces E t labeled by coordinate time 
t {— x° in our lab frame coordinate basis). The summary 
below serves to establish notation and spells out only the 
results we need here. Pedagogical introductions include 
Refs. dl-im. 

Generic metric components in the 3+1 formulation are 
found from consideration of a 'thin sandwich' of space- 
time bounded by two spacelike slices E t and Y, t +dt- In 
particular we consider the spacetime interval ds between 
two points: 2+ in E t , and 2+ + da+ in T> t +dt- With proper 
time interval dr orthogonal to E t , and proper length 
interval dl tangent to E t , ds is given by a Lorentzian 
version of the Pythagorean theorem (note the signature 
-+++): 



ds 2 = -dr 2 



(63) 



Denote the orthogonal proper time at a; M between E t and 



dr = a dt, 



(64) 



and call a the lapse function. In considering the proper 
length dt between x p in Et and the orthogonal projection 
of 2+ + G?a+ in Ti t +dt onto E t , we note that the curves 
x l — constant traced out by the spatial coordinates need 
not be orthogonal to E t . That is, the spatial coordinates 
may be moving as seen by an observer at rest in E t , such 
that they shift by a coordinate distance P % dt between E t 
and E t+( 2 t . Allowing for such a shift vector P % (which is 
tangent to E t ), and letting 7^ denote the three-metric 
within E t , we have 



dl 2 



jij (dx l + P l dt) (dx j + dt) 



(65) 



Comparing Eqs. (f6"3"|) - (f6"5)) with the line element ds 2 — 
g pv dx^ dx v ', we read off the metric components 



PkP k 



The inverse metric is 



9 



ft- 1 J 



-l/a 2 

fc/a 2 



7 



P 3 /a 2 
■ ' - P 1 P 3 ja 2 



(66) 



(67) 



The three-metric 7^ and its inverse 7 y lower and raise 
the indices of three- vectors within (tangent to) the space- 
like slice, as in /3j = 7y J $ J ' and ft 1 = ^ Pj- Finally, 



(68) 



expresses the determinant g of the four-metric in terms 
of the lapse function and the determinant 7 of the three- 
metric. 
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Solution of the Einstein equations for the met- 
ric components — nonlinear partial differential equations, 
second order in space and time — constitutes knowledge 
of spacetime. In the 3+1 approach, solution of the Ein- 
stein equations is transformed into a Cauchy problem: 
specify initial data (satisfying certain constraints from 
the Einstein equations) on an initial spacelike slice; and 
with coordinate freedom fixed and spatial boundary con- 
ditions specified, evolve the geometry of the spacelike 
slices forward in time. Of the ten degrees of freedom 
associated with the metric components (the number of 
independent elements in a 4 x 4 symmetric matrix), four 
correspond to the freedom in general relativity to choose 
any spacetime coordinates whatsoever, leaving six phys- 
ical degrees of freedom to be determined. The standard 
way to think about the coordinate freedom in the 3+1 
approach is to regard the lapse function a and shift vec- 
tor ft 1 as freely specifiable functions in time and space, 
associated respectively with the choice of time coordinate 
(i.e. the foliation or spacetime slicing) and the choice of 
spatial coordinates (in particular the motion of these co- 
ordinates as seen by an observer at rest in a slice). In 
order to facilitate practical solution by forward evolution 
in time, the second-order-in-time system is transformed 
to double the number of equations, first order in time, for 
double the number of dynamical variables. In particular, 
the above phrase "evolve the geometry of the spacelike 
slices forward in time" corresponds to evolution of (a) 
the six independent components of the three-metric 7^- 
governing the geometry within a slice, and (b) the six 
independent components of the extrinsic curvature Kij 
(another symmetric tensor tangent to the spacelike slice) 
that describes the warp of the spacelike slices as embed- 
ded in spacetime. Here we do not show these evolution 
equations; in what follows, we simply regard 7^ and 
as given, for instance as having been obtained by numer- 
ical solution (often of even further transformed systems, 
as for instance in BSSN and related approaches); see e.g. 

Refs. mum. 

B. Four- velocity of Eulerian observers and the 
spacetime divergence 

In dealing with various forms of stress energy and the 
equations that govern them in the 3+1 context, two help- 
ful tensors are the unit normal n M and the orthogonal 
projector 7^. The unit normal to a spacelike slice at 
a given point can be regarded as the four-velocity of an 
Eulerian observer, i.e. one at rest in the lab frame. In 
the lab frame coordinate basis its components are 

(n") = Q./a,-0 i /af, (69) 
(n M ) = (-a, 0,0,0). (70) 

Note that indeed n^n^ = g fll ,n p n L ' = — 1 as expected of 
a unit vector. The orthogonal projector is 

7*w = g^u + n^n v . (71) 



From Eqs. (|66| and ([70)) it follows that the spatial part 
of 7 M „ equals the three- metric 7^ , motivating use of the 
same base symbol. While contraction of an arbitrary 
vector with n p yields the portion orthogonal to a space- 
like slice, contraction with 7^ = g^ v + n^n v yields the 
portion tangent to the spacelike slice. Indeed a trivial 
calculation confirms that r ) il _ v n v = 0. 

The unit normal and orthogonal projector can be used 
to decompose a stress-energy tensor T^ v . The 'Eulerian 
projections' 

E = n^T^, (72) 
F» = -n v ^ p T v P, (73) 
= 7 M p l v ' a T pa (74) 

are respectively the energy density, momentum density 
(or energy flux) , and stress measured by an Eulerian ob- 
server. The momentum density and stress are spacelike, 
i.e. tangent to the spacelike slice: 

= 0, (75) 
n^S plJ = rivS^ = 0. (76) 

In terms of Eqs. ([72 )1 - ([73] ) . a stress-energy tensor can be 
decomposed as 

T"" = En p n v + F p n v + F u n" + 5"". (77) 

We call this the 'Eulerian decomposition' of T' 1 ". Em- 
phasizing the spacelike character of F p and S 1 *" (see 
Eqs. ([75"| and ([76]) ). it can be re-expressed 

T pv = E n^n v + F*'f i n u + X + S h V \-f ) (78) 

in the lab frame coordinate basis. 

Turning from a generic stress-energy tensor to the 
neutrino radiation in particular, we similarly define the 
Eulerian projections and Eulerian decomposition of the 
monochromatic neutrino stress-energy 

= L^L^T^, (79) 

whose comoving frame components T pv were given in 
Eqs. ([35]) and (|36|) . The coordinate transformation L^p, 
from the orthonormal comoving frame to the lab frame 
coordinate basis was discussed in Sec. Ill Al The Eulerian 
projections 

£ = n^T^, (80) 

^ = -n t ff P r' p , (si) 

= l\l v a T p ° (82) 

are respectively the monochromatic neutrino energy den- 
sity, momentum density, and stress as measured by an 
Eulerian observer. The momentum density and stress 
are spacelike, 

= 0, (83) 
n^ v = n v S» v = 0, (84) 
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and 



(85) 



is the Eulerian decomposition of the monochromatic 
stress-energy. 

The results in Appendix [A] for the four components of 
the spacetime divergence V M T'" / of a stress-energy ten- 
sor can be adapted immediately to the spacetime diver- 
gence (St) u of the monochromatic particle stress-energy 
in Eq. (|4ip. In particular we make the formal replace- 
ments 



E^£, 
F p -> P, 



(86) 
(87) 
(88) 
(89) 



From Eq. (|A14I) . the projection of the spacetime di- 
vergence orthogonal to the spacelike slice — which con- 
tributes to the neutrino energy equation — is 



n v (Sq 



d(D T>n ) , d (F r ,, 



Of 



dx l 



where 



(Ft,„) 1 



= SY(aJ*-P i £) 
P da 
a dx l 



-«V7 



■S %3 Ki 



(90) 

(91) 
(92) 

(93) 



are respectively the 'conserved' energy density, energy 
flux, and gravitational power associated with the neutri- 
nos. From Eq. (|A31|i . the projection of the spacetime 
divergence orthogonal to the spacelike slice — which con- 
tributes to the neutrino momentum equation — is 



lju (St)" = 
where 

(Ft, 7 )V 
(Gt )7 L 



1 



a(D Ti7 ). d(F Tn y 



dt 



dx l 



FT* 



(94) 

(95) 
(96) 



-(£ da 



2 dxi 



(97) 

are respectively the 'conserved' momentum density, mo- 
mentum flux, and gravitational force associated with the 
neutrinos. Note that £ (x p , e), T p (x p , e), and S» v (x p , e) 
are functions of the lab frame spacetime coordinates 
x p and the neutrino energy e reckoned in a comoving 
frame. Their relations to the primitive variables J (x p , e) 
and ^(x^e) follow from Eqs. P? ]) -(P |) . ([791. ©, and 
(f36|k but see also the following two subsections and Ap- 
pendix [Bj The projections of the spacetime divergence 
presented here are in accord with the corresponding 
terms in Eqs. (3.37) and (3.38) of Shibata et al. @. 



C. Four- velocity of Lagrangian observers and the 
momentum space divergence 

Before turning to the momentum space divergence, 
it will be helpful to begin rewriting in covariant form 
some of the expressions we have given involving comov- 
ing frame components, by writing them in terms of the 
four-velocity u M of Lagrangian observers (i.e. those at 
rest with respect to the fluid). In an orthonormal co- 
moving frame we have 



(1,0,0,0)' , 
(-1,0,0,0). 



Thus for example 



(98) 
(99) 



(100) 



is the neutrino energy measured by a Lagrangian ob- 
server, expressed in terms of lab frame coordinate basis 
components = upL p fl and p p = U l pp p (the coordi- 
nate transformation L M p, from the orthonormal comoving 
frame to the lab frame coordinate basis, and its inverse 
L p M , were discussed in Sec. Ill Al) . 

Turning to the comoving frame angular moments J , 
W ', and /C y " defined in Eqs. (I2"71) - (j3"0")) . we can define co- 
variant versions by extending t to a unit four-vector 
£ p satisfying £^l p = 1. This vector is spacelike in 
the comoving frame, with comoving frame components 
= (0,^W 3 ) T , and is thus orthogonal to it**, i.e. 



= 



(101) 



in any basis. Thus e£^ is a covariant representation 
of the three-momentum measured by a Lagrangian ob- 
server. Covariant representations of the comoving frame 
angular moments are 



J = e fdfl, 



W 



KT = 



(102) 
(103) 
(104) 
(105) 



where the integration is still performed with respect to 
the comoving frame solid angle. 

While in Eqs. (l3~Tj) - (l34l) we nominally defined closure of 
the system in terms of the comoving frame components 
1Q3 and C lJ of the second and third angular moments, 
in practice we may be able to obtain the lab frame co- 
ordinate basis components K. pv and C pvp without ever 
explicitly transforming any of the moments to and/or 
from the comoving frame using the transformations L p ^ 
and L p p discussed in Sec. Ill Al (Applying the closures 
without any reference to comoving frame components is 
desirable because it would be a hassle to have to work 
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explicitly with the transformations LP ^ and LP^. These 
will be quite complicated in the general multidimensional 
case, because neither the metric g^u, nor therefore the 
tetrad e M ^, nor either the boost A^, will be diagonal. 
Moreover, it would be necessary to apply the tetrad to 
the coordinate basis velocity variables obtained with the 
hydrodynamics solver in order to get the velocity parame- 
ters appearing in the boost. Working only with lab frame 
coordinate basis components of the comoving frame an- 
gular moments avoids these complications.) Following 
Ref. [44[ , we note that if K.p u is regarded as a function of 
the zeroth and first moments J and HP, the most general 
symmetric expression that (a) is spacelike relative to u p , 
and (b) satisfies the trace condition K,p M = J (obvious 
from Eqs. (fT02]) and (|T04| ). is of the form 

- -Jh^ + a(J, H) (WU V - ^H 2 h p A , (106) 



where H = y/W^H^, and 

hpw = g^u + u^Uy (107) 

is the orthogonal projector relative to u p . We extend 
this thinking to the third moment, finding the unique 
expression 

C^p = - (Whvp + n p h pu + n v h p p) 
5 



-h{J,U) 



WWW 



1 



■R 2 (HPh" p + H p hP" + % v h pli ) 



(108) 



satisfying the trace conditions (obvious from Eqs. (|103[) 
and (TUg)) 



rp f 



n p , cp\ = h v , c p \ = hp. 



(109) 



The point is that when KT" and C p " p are obtained only 
from knowledge of J and H p , in principle the only free- 
dom in the closures are the scalar functions a{J ,%) and 
b{S,%), in which the magnitude % = ^YW^W* can be 
evaluated in terms of lab frame coordinate basis compo- 
nents. 

Consider next the monochromatic stress energy T pv ■ 
As can be confirmed in the comoving frame, the neutrino 
momentum pp can be covariantly decomposed as 



e K + t p ) , 



(110) 



i.e. into portions tangent and orthogonal to the La- 
grangian observer four- velocity u p . Using this in the 
monochromatic stress energy 

T MV = ~ J pfpffdn (111) 

(see Eqs. (|35|) and ([79])). one can see that the 'Lagrangian 
projections' 

J = u^u„T pv , (112) 
W = -u„h%T y P, (113) 



do in fact yield Eqs. (|102l) - (|104[) . (Recall that h^ v , given 
by Eq. (|107[) . is the orthogonal projector relative to u p .) 
Thus the moments J , 7i p , and K.p v are respectively the 
monochromatic neutrino energy density, momentum den- 
sity, and stress measured by a Lagrangian observer. Note 
that W and K,p v are spacelike in the comoving frame: 



u^W = 0, 
uJO tu = u v K,p v = 0. 



(115) 
(116) 



Therefore 1~p v can be written as the 'Lagrangian decom- 
position' 



TP" = JuP 



(117) 



which provides an alternative to the Eulerian decompo- 
sition of Eq. (|85[) in terms of the energy density, mo- 
mentum density and stress measured by an Eulerian ob- 
server. 

Eqs. (|85p and (| 1 1 T[) can be used to write the monochro- 
matic energy density momentum density, and stress mea- 
sured by an Eulerian observer (£, TP, and SP") in terms 
of their counterparts measured by a Lagrangian observer 
{J, HP, and K.p v ), and vice- versa. In one direction, use 
Eqs. f8"0 ]) -([8"2" |l and substitute Eq. (fTTT)) on the right-hand 
side. In the other direction, use Eqs. (|112l) - (|114[) and sub- 
stitute Eq. (|85p on the right-hand side. We will say more 
about this in the next subsection and in Appendix [Bj 

Next we turn to the third momentum moment U p P v , 
given by 



UP vp = - / ,/■',/>''/,/<> 



(see Eq. ([37])). Note that 



Up U p P v = ~u p UP pv = -u p UP"P = (TP V 



(118) 



(119) 



is a covariant version of Eq. (|3"5|) . Using Eq. (|110[) in 
Eq. (fTT8|) and comparing with Eqs. ([Hj2 >(fT05)) , we find 



(J = -u^u v u p UP vp , 
e%P = u v u p hP a U vpo , 
cJCP" = ~ Up hP a h" ' K U paK , 
eC pup = hP a h u K h p X U™ X 



(120) 
(121) 
(122) 
(123) 



for the Lagrangian projections of U p,vp '. The associated 
Lagrangian decomposition is 

W v P = e ( Ju p u v u p + y_P-u v u p + H v uPu p + U p uPu v 

+K,p v uP + KPPu v + K pv uP + C,p vp ) . (124) 



K pv = hP p h\T 



■pa 



(114) 



Eqs. (|119[) and (|124j) are special to contraction and de- 
composition respectively with respect to the Lagrangian 
observer four- velocity u p . That is, the identity (up to a 
factor e) of the Lagrangian projections J, H^, and KT" 
of T M " with the Lagrangian projections of UP vp (except 
of course for the irreducible C,P vp ) holds due to our choice 
to measure energies and define angular moments in the 
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comoving frame. We can define (and in fact will use) the 
Eulerian projections oiU p " p , 

Z = -n^ripU^P, (125) 

y ,t = nyn p 'f a U vpa , (126) 

X^ v = -n pl ^ 1 \U paK , (127) 

W^ = r a7 \Y x ir K \ (128) 

and construct the associated Eulerian decomposition of 

U pvp = Zn^n v n p + yP n u n p + y v n P n p + y p n p n v 

+X pv n p + X w n u + X pv n p + W^ p . (129) 

The relationships between the Eulerian projections Z, 
y p , X pu , and W pvp of U^f and the Eulerian projec- 
tions £, J- p , and S pu of are not as simple as the 
relationships between the coefficients of Eqs. (|117j) and 
(|124|) . Nevertheless, useful relationships between the Eu- 
lerian projections of U p,vp and T^" do exist and will be 
presented in the next subsection. 

In a different vein, another comoving frame expression 
that can be written in covariant form thanks to Eq. (|99[) 
is 

= — T^^Uq — —^fipUp' = Vpiip — dj^uc, 

= V A ^, (130) 

a covariant expression for the connection coefficients ap- 
pearing in Eq. (|55]l. 

Having obtained these expressions rewritten covari- 
antly in terms of the Lagrangian observer four-velocity 
u p , we are ready to consider the momentum space diver- 
gence. Using Eq. (|130[) , the angle-integrated momentum 
space divergence of Eq. (|45l) can be written 

(M T ) P = ~ (-e 2 IfpUMVfiUe) , (131) 

or, taking advantage of the covariant nature of this ex- 
pression, 

(M T r = ^(-e 2 ^V^), (132) 

where U pvp is given by Eqs. (fT24|) or (fT29|) . Projecting 
orthogonal and tangent to the spacelike slice, we have 

- n p (M T ) P = ~ (e 2 npU^VyVv) , (133) 

Hp (Mt) p = ~ (-e 2 7j >W^V^) . (134) 

Up to multiplicative factors of e in defining the moments, 
the projections of the momentum space divergence (i.e. 
energy derivative) presented here are in accord with the 
corres pon ding terms in Eqs. (3.37) and (3.38) of Shibata 
et al. |20|. 



D. Three- velocity of Lagrangian observers and the 
relation between the lab and comoving frames 

In Section |n] we obtained expressions in terms of lab 
frame coordinate basis spacetime position components x p 
and comoving frame orthonormal basis momentum com- 
ponents p l through application of coordinate transfor- 
mations JJ 1 p,, but in the context of the 3+1 formulation 
there is a more fruitful way of expressing the relation- 
ship between the lab and comoving frames. In particular 
the Lagrangian observer four-velocity u p (i.e. the four- 
velocity of the fluid) can be 'Eulerian decomposed' into 
parts orthogonal and tangent to the spacelike slice, that 
is to say, parts tangent and orthogonal to the Eulerian 
observer four- velocity n M : 

u p = A(n p + v^) . (135) 

The orthogonality requirement on v^, 

n, L v p = 0, (136) 

implies (see Eq. ([70)) ) that v p is spacelike and has com- 
ponents 

K) = (0y) T (137) 

in the lab frame coordinate basis. The interpretation of 
v p as the three- velocity of a Lagrangian observer as mea- 
sured by an Eulerian observer is confirmed by squaring 
Eq. (|135[) to find the expected Lorentz factor 

A = (1 - v p v p y 1/2 = (1 - v'vi) ~ V2 (138) 

for a boost between frames with relative three-velocity 
v p . (Recall that u^u p = n^n^ = — 1. The scalar Lorentz 
factor A lacks indices and will not be confused with the 
Lorentz boost A^.) 

In the previous subsection we discussed finding the 
monochromatic energy density, momentum density, and 
stress measured by an Eulerian observer [£, J- p , and 
<S M ") in terms of their counterparts measured by a La- 
grangian observer [J, H 11 , and /C A " y ). In substituting 
the Lagrangian decomposition of Eq. (|117l) in Eqs. (|80|) - 
(JH3>, the factors of n p and 7^ „ can be expressed in 
terms of u p and v p via Eq. f|135|) . Using also Eqs. (|115[) 
and (|116p . this provides an alternate route to the rela- 
tions obtained more tediously from the transformation 
rj-tut = LtuVpTi* and Eq. (JMD ■ The explicit results 
are analogous to those obtained .271] via Lorentz trans- 
formations in special relativity. They are not particu- 
larly illuminating, but for completeness we exhibit them 
in Appendix IBl Relations between £, F p 1 S pi/ and J, 
H p , K,^ will of course be needed in numerical work, but 
it may be best to perform the contractions in Eqs. (|80p- 
([82]) numerically rather than code the tedious analytic 
expressions in Appendix |B| 

We also use Eqs. (|135[) and (|136|) in fulfilling our 
promise, made in the previous subsection, to relate the 
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Eulerian projections of U^ p in Eq. ([129]) to those of 
in Eq. (j8"5)l . These relations are obtained by plugging 
Eq. (|129j) into Eq. (|119|) and comparing the results with 
Eq. (1551) for the coefficients of outer products of two, one, 
and zero copies of n p . The results are 

A(Z-tv^) = e£, (139) 
A (y" = cF", (140) 
A (A 1 ^ - v p W^ p ) = eS M ". (141) 

These can be 'unraveled' in reverse order to give 

A^"= £ 5^ + KvpW^P, (142) 
A = e J* + v„ (e 5"" + Kv p W pup ) , (143) 
A Z = e £ + [eJ* + v v (e 5"" + A^W"")] . 

(144) 

We discussed finding the Eulerian projections £ , J- 71 , and 
5' tl/ in terms of the Lagrangian J, W 1 , and TC^ V (i.e. 
the comoving frame angular moments) in the previous 
paragraph; see also Appendix [Bj There remains the third 
moment W A " yp , which can be found in terms of J , W p , 
K,^ , and C^ vp through the projection 



(TK.X 



(145) 



using the Lagrangian decomposition of Eq. (|124[) on the 
right-hand side. Indeed this same approach can be used 
in lieu of Eqs. (|T42l) - (fT44)) to directly obtain Z, JW*, and 
X 1 *" in terms of J, W, W , and L pup as well. Again 
we reserve explicit expressions for Appendix IB] 

These Eulerian projections come into play in mak- 
ing the momentum space divergence more explicit, for 
we use the Eulerian decomposition of U pupj given by 
Eq. (|129l) and the Eulerian decomposition of u v given 
by Eq. (|135l) in expanding the expression U pvp V ^u^ ap- 
pearing in Eqs. (| 1331) and (|134l) . At first glance the use 
of these Eulerian decompositions may seem to compli- 
cate things. With our choice to measure neutrino mo- 
mentum components and define angular moments in the 
comoving frame, the Lagrangian decomposition of U pvpj 
in Eq. (I124[) is simpler in the sense described in Section 
IIII C[ and indeed u v is the four-velocity of Lagrangian 
observers; so why not stay with these Lagrangian expres- 
sions? The problem is that we cannot stay in 'Lagrangian 
world' — in the comoving frame — altogether, much as we 
might like to, because the covariant derivative V^u,, is 
with respect to the lab frame coordinate basis. We have 
swept the connection coefficients in Eq. (|4"5"|) temporarily 
under the rug via Eq. (|130p . but lab frame coordinate 
basis connection coefficients still lurk in the covariant 
derivative W ^u u . 

Given the unavoidable necessity of facing the relation 
between the lab and comoving frames in one way or an- 
other, there are significant advantages to consistent use of 
Eulerian decompositions, in which the relation between 
frames is focused more in the three- velocity of a La- 
grangian observer as measured by an Eulerian observer, 
than in the coordinate transformations L p p, of Eq. (|7J). 



If we stay with Eq. (|45|) . we face the unpleasant 
prospect of evaluating Eq. (TTT]) for the transformed con- 
nection coefficients. Even rewritten as Eq. (II 32[) . we face 
lab frame coordinate basis connection coefficients when 
using the Lagrangian decomposition of U pvp and leaving 
the Lagrangian observer four-velocity u v as is. But as 
derivations in Appendices [S] and [C] show, the fact that 
Eulerian decompositions are most natural in the 3+1 ap- 
proach allows us to almost completely avoid explicit en- 
counters with connection coefficients. 

Moreover, consistent use of the Eulerian perspective — 
both in projecting out the portions of the phase space 
divergence orthogonal and tangent to the spacelike slice, 
and in Eulerian decompositions of U pvp and u v — also 
turns out to preclude any appearance of time derivatives 
of metric functions, even in intermediate steps. Time 
derivatives of the lapse function a and shift vector (3 l 
would be particularly inconvenient in numerical work, as 
these do not normally have evolution equations associ- 
ated with them. (Unfortunately, we shall see that time 
derivatives of the Lorentz factor A and three-velocity v l 
remain; these are something of a nuisance, but at least in 
principle they could be written in terms of spatial deriva- 
tives via hydrodynamics evolution equations.) 

Finally, the Eulerian decomposition of U pvp " is more 
readily tied to the Eulerian decomposition of T 1 ^, i.e. to 
the energy density £ , momentum density T p , and stress 
S pu measured by an Eulerian observer. This is advanta- 
geous in relating four-momentum conservation to lepton 
number conservation, for cancellations must occur be- 
tween the spacetime and momentum space divergences, 
and it is £ , J- p , and that appear in the spacetime 
divergence in Eqs. (|9"0)) and 

Details of the calculation of U pup V ^u v , using the Eu- 
lerian decompositions of U pvp and V M Uj,, are given in 
Appendix [C] Using those results for the projection or- 
thogonal to the spacelike slice, we have 



n u (M T ) v = 



1 1 d 



[e 2 (Rt> + Or,*)] - (146) 



where 



-X kl v m 



a dx l dx l 



dx r < 



X kl K kl 



dx^ ~^ dx l 



-y k n» d -^-x/^ 



dx^ 



dx % 



(147) 



(148) 



arise from changes in the neutrino energy as measured in 
the comoving frame due to gravitational redshift (Rt.m) 
and the acceleration of the observer riding along with the 
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fluid (Ot.h)- Similarly, 

1 1 d 



l3» (MtT 
with 

(Rt, 7 ), = 

(0t, 7 ), = 



{e 2 [(R T , 7 ) . + (0 T)7 )J } , 

(149) 



e 2 <9e 
0>' - X/) da 



a 



X n W 
ftr* ife ftr* 



9 (Av k ) 



dA 



(150) 



(151) 



for the projection of the momentum space divergence tan- 
gent to the spacelike slice. That the structures of Rt. m 
and Ot.u parallel those of (Rt, 7 )^ and (0t )7 )^. is sim- 



ply a reflection of the parallel structure of Eqs. (|C3I) and 
(|C4[) . As given here they are expressed in terms of the 
Eulerian projections Z, y», X» u , and W^ p oiU pv ^ (see 
Eq. (|129jl ). In numerical work these could be further ex- 
pressed in terms of the Eulerian projections of T" M (see 
Eq. (18510 . i.e. the energy density £ , momentum density 
J 7 ' 1 , and stress S pi/ measured by an Eulerian observer, 
via Eqs. (|142l) - (|144|) . They could also be expressed di- 
rectly in terms of the Lagrangian projections, as shown 
in Appendix [Bl 



E. Four-momentum and lepton number exchange 

The relationship of lepton number and energy ex- 
change at a high level is readily seen. Revisiting Section 
III Dl in terms of covariant expressions involving the La- 
grangian observer four- velocity w M (see Section IIII C|) , it 
is easy to see that Eq. (j5"Tj) becomes 



A/ 7 " = MuT^y 

e 



that Eq. (|54j) becomes 

1 d 



(152) 



(153) 



and that Eq. (160)) becomes 



1 



Sat + M w = -~ Uu [(S T y + (M T ) v ]. (154) 



The first term on the right-hand side is 

^ (155) 



-u v (St)" = --u^V^T^ 



( 1 \ T M " 
V M [--UuT^J + — V^u„(156) 

T-jUI/ 

S w + V„Uu, (157) 



thanks to Eq. f| 152|) . The second term on the right-hand 
side of Eq. (|154|) is 



-u p (M T y 
e 



1 

--u p 
e 

l_d_ 

e 2 de 



e 2 T^V^Uv 

8(1 
e 



^ v dt 



>N 



(158) 

(159) 
(160) 



by virtue of Eqs. (1TT§| and (1T53")) . The sum of Eqs. (fT5?)) 
and (|160p gives Eq. (|154|) as required. 

This consistency between lepton number and energy 
exchange applies not only at this high level, but also to 
the detailed form of the moment equations as we have 
most expressly written them, which has implications for 
their discretization. In terms of the decomposition of 
in Eq. (I135[) . we have 



5JV 



Mat = -- (n, 



+ v u )[(S T ) v + (MtH (161) 



for Eq. (|154[) . This is, naturally, closely related to our 
projections of the spacetime and momentum space di- 
vergences orthogonal and tangent to the spacelike slice. 

We consider first the spacetime divergence on the right- 
hand side of Eq. (fTBTT) . Using Eqs. fl90j)-([92]) and (j69l . 
we have 



^-n„ (St)" 





[A(D T ,„)" 


d 


A(F T ,„) r 


\dt 


e 


dx l 


e 



[A(G Tl n) 



where 



-T,n 



(162) 



(163) 



are the 'extra' terms that arise from pulling the factor 
A inside the time and space derivatives. Thus the dis- 
cretized form of Et,h will be dictated by the discretiza- 
tion chosen for the first two terms of Eq. (|90l) . Similarly, 
using Eqs. ([M])-(|M]) and we have 



7jv (S-2 



a^' (d Ti7 ) : 


d 






dx l 


e 



Av j (G T , 7 ), + E T/ 



1 



where 

ET i7 = Cty/j 



/i d(Av 3 ) gi d(Avi) 



dx p 



(164) 



(165) 
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are the 'extra' terms that arise from pulling the factor 
Kv 3 inside the time and space derivatives. Thus the dis- 
cretized form of Et, 7 will be dictated by the discretiza- 
tion chosen for the first two terms of Eq. Adding 
Eqs. (|T62l) and (H"64| gives 



( n v V^v) (St)" 
A 



Sjv - S — 
1 



(G*r>) - v 3 (G T , 7 )j 



(166) 



for the contribution of the spacetime divergence to the 
right-hand side of Eq. (|161l) . 

Turning to the momentum space divergence, from 
Eqs. (I146P and (|149l) we have 

- ^ ( n v + u^v) (M T )" 



ijv 



1__5 

i 2 ^ 

+Ae 2 



{Ae 2 [(R T ,„) - v 3 (R Tn ) J 
(Or, n ) - ^ (0 T , 7 ) 1 } J- (167) 



for the contribution of the momentum space divergence 
to the right-hand side of Eq. (|161[) . Note that the long 
term following Mjy results from pulling 1/e through the 
energy derivative; therefore its discretized form is dic- 
tated by the discretization chosen for Eqs. (jl46[) and 

CHS) . 

The sum of Eqs. {166)) and ([T6T|t equals Eq. (fTBTjl . 

as required, because (a) the gravitational redshift terms 
from the momentum space divergence cancel the grav- 
itational force and power terms from the spacetime di- 
vergence, and (b) the observer corrections from the mo- 
mentum space divergence cancel the 'extra' terms from 
pulling A and v k through the time and space derivatives. 
That is, (a) 



0=* 
e 

and (b) 



(Gt>) - v 3 (G T , 7 ) 



A 

72 



(168) 







\ (Et,, 



E T, 7 ) + 4 



(O r ,n) - v j (0r l7 )j 



(169) 



These cancellations emerge in a surprisingly tractable 
way: the gravitational redshift and observer correction 
terms from the normal and tangent projections of the 
momentum space divergence combine in just the right 
way to make use of Eqs. (|139j) - (|141|) . which relate the 
Eulerian projections Z, y, ', and of W pl ^ 

(see Eq. (|129[1 ) to the Eulerian projections of T vti (see 
Eq. (|85jl). i.e. the energy density £ , momentum den- 
sity J r/i , and stress iS Ml/ measured by an Eulerian ob- 
server. In particular, combining Eqs. (|147|) and (|150[) 



using Eqs. (fT39|) - (fT4T|) yields 



-S fc V 



dx % 
- S kl Kki 



d(3 k 
dx % 

(170) 



Plugging this and Eqs. flgSJ) and (JST) into Eq. ([I55| . the 
term- by-term cancellations are apparent. Similarly, com- 
bining Eqs. (fT48"l) and (|T5T|) using Eqs. (fH9>(fT4Tj) yields 



(O r ,n) ~ v 3 (0 Tn ) 3 



A 



«V7 



gA 



P 



dk 

dx l 



-Fun 1 



,d(Ai 



dx^ 



S k 



(171) 



Plugging this and Eqs. (fToU)) and (fTo5)) into Eq. (ITM| . 
once again the term-by-term cancellations are apparent. 

This analytic demonstration of the consistency of our 
conservative four-momentum moment equations with a 
conservative number moment equation ideally should be 
repeated in the discretized case in order to discover dis- 
cretizations that are faithful to this consistency. 



IV. CONCLUSION 

We now assemble the expressions obtained in 
Sees. MI BllII Dl into conservative 3+1 general relativistic 
variable Eddington tensor radiation transport equations. 
These four equations are conservation laws for the energy 
and momentum carried by the neutrino radiation: 

d(D Tn ) d(F Tn Y 1 d r 2 ,„ n Nl 

= G T ,« + C T ,n, (172) 
d(D TtJ ) j c)(F T , 7 )^ i o 



Of 



dx l 



de 



{e 2 [(R T)7 ) . + (0 T)7 )J } 



(Gt, 7 L + (C T , 7 ), 



(173) 



The global 'lab frame' spacetime coordinates t and x l 
are those associated with the 3+1 formulation of general 
relativity, in which the line element and metric compo- 
nents g^ v are given by Eqs. (jB"3")) -([6"5 j) . Equations (|172[) 
and (|173p . expressing energy and momentum conserva- 
tion respectively, come from the projections of the phase 
space divergence of the monochromatic neutrino stress 
energy orthogonal and tangent to the spacelike slice. The 
projection orthogonal to the spacelike slice is via con- 
traction with n^, the unit normal to the spacelike slice, 
which is also the four- velocity of Eulerian observers (see 
Eqs. (|6"9l and (1701) ). The projection tangent to the space- 
like slice is via contraction with the orthogonal projector 
liiv = 9 'nv + n^n u . The momentum space coordinate in 
these angle-integrated moment equations is the energy 
e = —u^p 11 measured by a Lagrangian observer, whose 
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four-velocity it M is that of the fluid. These coordinate 
choices allow particle/fluid interactions to be evaluated in 
the comoving frame in the context of Eulerian grid-based 
approaches to multidimensional spatial dependence. 

A rough analogy with conservative formulations of hy- 
drodynamics is evident. The 'conserved' monochromatic 
energy and momentum densities n and (Dy~) . are 
given by Eqs. (|91l) and (|95[) . The energy and momentum 
fluxes (F*r,n) 1 and (Ft !T )'. are given by Eqs. (li?2"j) and 



(|96p . These are expressed in terms of the 'Eulerian pro- 
jections' (see Eq. (|85jl ) of the monochromatic neutrino 
stress energy T Miy — the second momentum angular mo- 
ment of the distribution function / (see Eq. (|111[) ). These 
Eulerian projections are the monochromatic energy den- 
sity £, the momentum density J-^, and the stress S^ v 
measured by an Eulerian observer (i.e. in the lab frame). 
These may be expressed (see Appendix [Bj) in terms of the 
'Lagrangian projections' of T^ v (see Eq. (|117|) ). which 
are the energy density J7, the momentum density W 1 , 
and the stress XT 1 " measured by a Lagrangian observer 
(i.e. in the comoving frame). The closure relations — 
that is, the Eddington tensors — are defined in terms of 
these Lagrangian projections. Roughly speaking, the Eu- 
lerian and Lagrangian projections are respectively like 
the 'conserved' and 'primitive' variables in hydrodynam- 
ics, while the Eddington tensor closure relation between 
Lagrangian projections is analogous to the equation of 
state relating primitive (comoving frame) hydrodynam- 
ics variables. 

As with self-gravitating hydrodynamics, the presence 
of gravitational source terms makes Eqs. (|172[) and (j 1 73[) 
more properly 'balance equations' rather than strict con- 
servation laws. The terms Gr, n and (Gt^)^, given by 
Eqs. (|93| and (|97|) . represent the energy and momentum 
exchange between the neutrinos and the gravitational 
field as embodied in the spacetime geometry. (If spher- 
ical or cylindrical spatial coordinates were used, these 
source terms also would include the fictitious forces as- 
sociated with these. For the purpose of discretizing the 
spatial flux, it would be desirable in this case to factor 
out the portion of the spatial metric determinant arising 
from these coordinate choices, such that it would appear 

(a) in the denominator outside the spatial derivative, and 

(b) in the numerator inside the derivative, separate from 
the flux.) 

As monochromatic (i.e. energy-dependent) radiation 
transport equations, Eqs. (|172[) and (|173l) contain terms 
beyond those present in self-gravitating hydrodynamics. 
First, the source terms 



(C T , 7 ). = aV77j>^ 



(174) 
(175) 



represent the energy and momentum exchange with the 
fluid, with the relationship of the source to the col- 
lision integral given in Eq. (|46|) . Second, the energy di- 
vergences on the left-hand sides of Eqs. (|172|) and (|173[) 
arise from changes in the neutrino energy as measured in 



the comoving frame due to gravitational redshift — Rt,u 
and (R.T n )j, see Eqs. (|147|) and (|150|) — and the acceler- 
ation of the observer riding along with the fluid, Ot,u 
and (Ot,j)j', see Eqs. (|148[) and (I15ip . These are given 
in terms of the Eulerian projections (see Eq. (|129[1 ) of 
the monochromatic third momentum angular moment 
U^p (see Eq. (fTTg)l ). The Eulerian projections of W vp 
can be expressed in terms of the Eulerian projections 
of T MI/ (Section IIII Dp . or more directly in terms of the 
Lagrangian projections (Appendix IB"]) . Finally, we note 
that the energy divergence is conservative with respect 
to integration over the differential energy volume e 2 de. 

While the conservative variable Eddington tensor mo- 
ment equations express four-momentum exchange be- 
tween the neutrinos and the fluid, we also have exam- 
ined in detail the relationship of these to the conservative 
number exchange equation (see Sections III Dl and IIIIEp . 
This can be done in a tractable way on a term-by-term 
basis thanks to our greater elaboration of the momen- 
tum space divergence than in previous work. Important 
conceptual features of our approach include (a) consis- 
tent use of what we call 'Eulerian decompositions' and 
'Eulerian projections,' which are natural to the 3+1 for- 
mulation; and relatedly, (b) a shift from conceptualizing 
the relationship between the lab and comoving frames 
from coordinate transformations L^p, to the (covariant) 
relative three-velocity connecting the four-velocities 
and u p of Eulerian and Lagrangian observers. Our 
approach is geometric, in conception if not notation, to 
an extent that allows us to obtain explicit results while 
almost completely avoiding encounters with connection 
coefficients. This understanding of the relationship be- 
tween conservative four-momentum exchange and conser- 
vative number exchange can guide the determination of 
discretizations of the variable Eddington tensor moment 
equations that facilitate simultaneous energy and lepton 
number conservation in numerical simulations, yielding 
greater confidence in simulation outcomes. 



Appendix A: Spacetime divergence 

In this appendix we compute the spacetime divergence 
of a stress-energy tensor with the 3+1 metric. The Eule- 
rian decomposition of a stress-energy tensor is given by 
Eq. I|77p or (|78p . In particular, we derive contributions to 
'conservative' evolution equations for the energy density 
E and momentum density Fi measured by an Eulerian 
observer (whose four-velocity is the unit normal n u to a 
spacelike slice). 

Some relations involving derivatives of the unit normal 
and the orthogonal projector will prove useful. 
The gradient of the unit normal is related to the extrinsic 
curvature and lapse function by (4l| 



V7 - -VT "m da 



(Al) 



Because n^V^n^ = V /J (n 1/ n 1/ )/2 = 0; and because K^ v 
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is tangent to the spacelike slice, i.e. spacelike in the lab 
frame coordinate basis (n M if M „ = n"Kp V — 0), the non- 
vanishing projections of this equation are 

1 da 



i dx u 
-Ku 



(A2) 
(A3) 



in the lab frame coordinate basis. Eq. ()A2j) relates the 
four-acceleration of an Eulerian observer to the gradient 
of the lapse function. Eq. (IA3[) relates the spatial part of 
the gradient of the unit normal to the extrinsic curvature, 
expressing the fact that the direction of the normal varies 
with the warp of the slice as embedded in spacetime. 
Another relation valid in the lab frame coordinate basis 



for vectors tangent to the spacelike slice (z° 



dn^ 



-Zi- 



This follows from writing 



(z M spacelike) 



dn^ 
'~dx^ 



zo 



dn° 
dx v 



goiz 



,dn° 



dri 1 
dri 



0) is 
(A4) 

(A5) 
(A6) 



dx v dx v 

and using goi = and Eq. (|69[) . Finally, and more 
straightforwardly, the gradient of the orthogonal projec- 
tor is 



(A7) 



thanks to the vanishing covariant derivative of the four- 
metric gp V . 

The projection of the spacetime divergence of a stress- 
energy tensor orthogonal to the spacelike slice contributes 
to an energy equation. Contracting the divergence with 
n v and taking it inside the derivative, we have 

- n v VpT» v = -V M (nj"") + T^V^. (A8) 
The first term on the right-hand side is 

- V M (n.TH = V M (En* + *V,) (A9) 



1 







'—g dx^ 
1 d 



/=g(Bn"+F i y'i)][A10) 



dt 
1 d 



I [^{aF i -l3 i E)] , (All) 



where we have used Eqs. ([M]). and (jTS]) . In substi- 

tuting Eq. (|78p into the second term on the right-hand 
side of Eq. (|A8[) . the first two terms vanish: 

(A12) 



because n^V^n,, 
terms give 



V t _ l (n v n v )/2 = 0. The remaining 



F i da 
a dx l 



(A13) 



thanks to Eqs. (|A"2j) and (|A3l . Putting Eqs. (|A11|) and 
(|AT3| together, 



Qfy^ dt 

_L d_ 

a \fl dx l 
F i da 
a dx l 



S ij K v 



(A14) 



is the portion of the spacetime divergence orthogonal to 
the spacelike slice. 

The projection of the spacetime divergence tangent to 
the spacelike slice, which contributes to a momentum 
equation, is a bit more involved. Contracting with the 
orthogonal projector and taking it inside the derivative, 

IjuV^ = V M (7;^"") - T^V^jV- (A15) 
The first term on the right-hand side is 
1 d 



-g dx^ 



The first term on the right-hand side is 
1 d 



(A16) 



-g dx^ 



( 



1 d r 



l—g dx^ 
1 d 



-g{F in * + &ff$[ (A17) 



dt 
1 d 



[yftia&i-pFi)] .(A18) 



a \/l dx 

The second term on the right-hand side of Eq. (|A16|) is 



-T" ml p^^-V m {Fp^ + S%). 
The first term on the right-hand side is 

- - F p {d 3 nP - V 3 nP) 

dft 



= -Fi 



dxi 



F l K 



(A19) 

(A20) 
(A21) 



where the first and second terms follow from Eqs. (|A4I 
and (|A3p respectively. The second term on the right 
hand side of Eq. (|A19|) is 



jp p 



-r£ g pk s 



■/ k 



(A22) 



= -?f?j + ^r-SV(A23) 

2 \ ox-? ax 4 ox' 1 / 

The last two terms in parentheses are antisymmetric in 
i, k and vanish upon contraction with the symmetric S , 
leaving 



■pp CM 

p 



S lk d llk 
2 dxi ■ 



(A24) 
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With Eqs. (|MT|) and (jA24l) . Eq. (jM9|) becomes 

-r^T- = -^g + F^,-f (A25) 
This, together with Eq. (IA18|) . yields 
V, ( 7j , T n = -^=|:(V7^) 



1 .9 



+ 



a^jFf dx l 



1 dx> + * 3 l 2 dxJ 



(A26) 



for Eq. (|A16j), th e first term of Eq. (|A15j) . The second 
term of Eq. (|A15|) is less complicated. Using Eq. (|A7|) . 

- T^V^ jv = -T^ (n v V„ nj + rijV^) (A27) 
= -T^nyVpnj, (A28) 



in which the second term has vanished due to Eq. (f70|) . 
This can be rewritten 

- T» v V^ jv = (En" + F^\) (A29) 
E da 



(A30) 



a dxi * K ^ 



where we successively have used Eq. ([75)1 and Eqs. (|A2[) 
and (|Al!j) . Adding Eqs. (|A26|) and (|A30I) and allowing 
for the symmetry of K^j , we have 



i a 



o 



»/7 <9x l 



[V7(«^--^-)] 



for the projection of the spacetime divergence tangent to 
the spacelike slice. 



Appendix B: Eulerian and Lagrangian projections 

In this appendix we express the Eulerian projections 
S, F p , and S py oiT^ (see Eqs. (JTTTJ and ([8§])). and Z, 
3^, A^", and of (see Eqs. (JTTBJ) and ([T^l) ). 

in terms the angular moments J, H p , /C M ", and 
(see Eqs. (|102[) - (|105[) ). which are also the Lagrangian 
projections of T pv and — up to a factor e — of U pvp (see 
Eqs. (I117[) and (|124[1 ). Because the Lagrangian projec- 
tions are spacelike in the comoving frame, we eliminate 
u p in favor of v p and (in free indices) n/ i . In the case of 
the neutrino energy density, flux, and stress, we also show 
how the relations are consistent with published results in 
special relativity obtained with Lorentz transformations 
27], rather than with projections as done here. 



We use four basic types of contractions, which follow 
from Eq. (fT3l)j) : 



UpU^ = A, 



-VpZ p , = 0) 



(Bl) 
(B2) 

7p X = Av p , (B3) 
7pm 2 " ^ fcw ~ «/>"/*) ^, = 0) (B4) 



where Z M is a stand-in for any spacelike index relative to 
the Lagrangian observer four- velocity u p (i.e. UpZ p = 0). 
The Eulerian projections of T pv are 



£ = n tl n v T IJ ' v ', 



Spa = Ip^lavT" 



(B5) 

(B6) 
(B7) 



Using Eqs. (fTT7) and ([BT ]) -(|B3 ]| on the right-hand sides, 
we find 



£ = A 2 J + 2AvpW + v^vJKT , (B8) 
T p = A 2 v p J + A (g pli - n p Vp) W 

+Av p VpW + (g P p - n pV p) v v K pv , (B9) 
S pcr = A 2 VpV a J + A (g P p - n p Vp) v a H p 
+A (g a p -n a Vp) v p W 
+ (5pm - n p Vp) (g„ v - n a v u ) K pv . (BIO) 
The Eulerian projections of U^ p are 

Z = -np L n v n p W l,p , (Bll) 



ya=la t ,n v n p U pup 1 

Xan = -la l i r lK V npU pl ' P , 



(B12) 
(B13) 
(B14) 



Using Eqs. (fi~24|) and (|Bl ]) -(|B4 )l on the right-hand sides, 
we find 

Z = e (A 3 J + 2,A 2 VpU p + ZAv^K,^ 

+v^v p £ p » p ,) (B15) 
y a = e [A 3 v a J + A 2 (g a p - n a Vp) W + 2A 2 V<r v^W 
+2A (gap - n a Vp) v u K pu + Av a VpVvK pv 
+ (ga^ - naVp) v v v p L pvp \ , (B16) 
Xan = e [A 3 v a v K J + A 2 (g a p - n a Vp) v K U p 

+A 2 (g K p ~ n K Vp) v a U p + A 2 v a v K VpH p 
+A (gap - riaVp) (g KU ~ n K v v ) K pv 
+A (gap - riaVp) v v K pv + A (g K p - n K Vp) v v K pv 
+ (gen - naVp) (g KV - n K v v ) v p Cy- vp \ , (B17) 
W ctkA = e [A 3 v a v K v x J + A 2 (g a p ~ n a Vp) v K v x H p 
+A 2 (g K p ~ n K Vp) v\VaW 
+A 2 (g\p - n x vp) VaV K H p 
+A (gap - riaVp) (g RV - n K v v ) v\K pv 
+A (g K p - n K Vp) (g Xv - n x v u ) v a K pv 
+A (gxp - n x vp) (g av - n a v u ) v K K pv 
+ (gap - n a Vp) (g KV - n K v v ) 

x (gx P - n x vp) C pvp ] . (B18) 
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We now prepare to compare the (raised index) spa- 
tial components of Eqs. (|B8I) - (|B10|) with the special rel- 
ativists results in Eqs. (182)-(184) of Ref. 0. In flat 
spacetime the unit normal of Eq. (|69p becomes 



K) = (1,0,0, of 



(B19) 



The tetrad is a Kronecker delta in flat spacetime 
Cartesian coordinates, so that the composite transfor- 
mation (x of Eq. ([7]) becomes 



(B20) 



An explicit expression for the Lorentz boost between the 
orthonormal lab frame and the comoving frame is 



A° 6 A° ? 
A*. A T | 



A 

AV ri 



(A-l) 



V l Vs 



where 

V 2 = V l V t + VV 2 + V S V§. 
It can be shown that 



(B21) 



(B22) 



A = 1 - V) 



-1/2 u s-1/2 /-. i 

= (1 - V^Vfj,) = ( 1 - V Vi 



-1/2 

(B23) 

is equal to the Lorentz boost we have been using. The 
quantities 



V 1 = V i , V 2 = Vi, 



V A = V, 



(B24) 



are not to be regarded as components of a four-vector, 
but simply as the three-velocity parameters appearing 
in the Lorentz boost, expressed in a manner consistent 
with our index conventions. In flat spacetime the spatial 
components of the (covariant) three-velocity four- vector 
u M are related to these boost velocity parameters by 



6\V\ 



(B25) 



a perhaps expected result that follows from Eqs. (|135l) . 
(|Bl9|) , (|B20|) , (|B2T]) . and flM]). However, we caution that 
the perhaps less-expected result for the lowered-index co- 
moving frame spatial components of is 

Vi = AV t . (B26) 

Below we use this to evaluate contractions of the form 



= KViZ\ (u^Z" = 0) 



(B27) 
(B28) 



i.e. contraction with an index that is spacelike relative 
to the Lagrangian observer four- velocity vf*. 

We are now ready to compare the (raised index) spa- 
tial components of Eqs. (IB8|) - (|B10I) with the special rel- 
ativistic results in Eqs. (182)-(184) of Ref. [27]], using 
Eqs. ([BT9]) , (jB20| . (|B2T|) . and (|B28|) . Equation dH5} can 
be expressed as 



£ = A 2 J + 2Av t W + v { v s K.v, 



(B29) 



£ = A 2 J + 2A 2 V{H i + A 2 ViV,fC i3 , 



(B30) 



which agrees with Eq. (182) of Ref. [27[. The spatial 
components of Eq. (|B9[) are 



F = A 2 v l J + AL\W + Av l v{W + L'tvjlC 3 , (B31) 



F = ^ |A 2 V\7 + A 
+A 2 V l V{H i + A 



W 



(B32) 



which agrees with Eq. (183) of Ref. [23|. The spatial 
components of Eq. (|B9|) are 



S ij = A 2 v i v j J + AL\v j W + \ I.' ,! ' It/ + /.',/. ' A". 

(B33) 



or 



5 IJ = S^j \ A 2 V l V 3 J + A [d-.Y'fi: + S^V'W) 



+2A 



r 2 



(B34) 



which agrees with Eq. (184) of Ref. |27j . 



Appendix C: Momentum space divergence 

In this appendix we compute, with the 3+1 metric, ex- 
pressions entering the momentum space divergence term 
of the phase space divergence of a monochromatic stress- 
energy tensor. In particular, we derive contributions to 
'conservative' evolution equations for the monochromatic 
energy density £ and momentum density Ti measured by 
an Eulerian observer (whose four-velocity is the unit nor- 
mal to a spacelike slice). 

We project the momentum space divergence into por- 
tions orthogonal and tangent to the spacelike slice, and 
find it helpful to also use Eulerian decompositions of the 
tensors that appear. Our starting points are Eqs. (|133[) 
and (|134[) for the portions of the momentum space diver- 
gence orthogonal and tangent to the spacelike slice: 

-n p (M T ) p = ^-^-(e 2 n p U^\7^u u ), (CI) 

Hp (Mt)' = ~ {-e 2 ljp U»^V»uJ) . (C2) 

The Eulerian decomposition of U pu ^ is given by 
Eq. (|129[) . The projections orthogonal and tangent to 



18 



the spacelike slice, which appear in the above two equa- 
tions, are 

~n p U plJfl = Zn v n p + y u n p + y p n u + X v », (C3) 
j jp U pvfi = y 3 nW + X 6 W + X/n u + W/ M (C4) 

We also use Eq. (fH5j) . 

= A (n^ + v p ) , (C5) 

i.e. the Eulerian decomposition of u v , in expressing 

There are four types of contractions appearing when 
either Eq. (|C3[) or (|C4[) is contracted with V p u„; we con- 
sider them in turn, beginning with 

nWV p u v = nW [(n„ + v v ) d^A + A V M n v +AV M i;„] . 

The first term gives 

nW (n v + v v ) d p A = -n"d„A (C7) 

because of the orthogonality relation n^v^ = 0. The 
second term vanishes because n p n^ = — 1 = constant 
implies 

= V M (n"n v ) = 2n"V„n„. (C8) 

In the third term we use 

n"V M «„ = -^n,, (C9) 

which follows from = V M (n v v u ). The third term be- 
comes 



n"^ (AV M «„) = - 



Av l da 



a dx % ' 

where we have used Eq. (IA2p . All together, 



(CIO) 



A^ 5a „ <9A 



for this type of contraction. 

Consider next a contraction of the form 



AWV^ = AW [n v d p A + AV M n„ + V M (Av„)] , 

(C12) 

where A v is spacelike (i.e. = 0, such that A = 

in the coordinate basis), and we break up the right-hand 
side differently than in the previous paragraph. The first 
term vanishes because A" is spacelike. The second term 
gives 



AW (AV,^) 



AA i da 
a dx l 



(C13) 



thanks again to Eq. (|A2[) . In the third term it will turn 
out best to raise the index on v v before turning the gra- 



dient back, in a sense, on n^: 

AW [V M (A«„)] = AW [V„ (Av k )] (C14) 
= AW[d^(Av k )+T h m ^Av m } (C15) 
= AWd„ (Av k ) 

+AA k v m (rj^n") (C16) 
= AWd„ (Av k ) 

+AA k v m (W m n k - d m n k ) (C17) 
= AWd„ (Av k ) - AA k v m K mk 

+AA k v m d m /3 k , (C18) 

where we have used Eqs. (|A3|) and (|A4[) in the last step. 
All together, 



AWV,u v 



AA* 8a . . % d$ k , 4 , 



^-7 -i- ~ AA v l K ik 

a ox ox 



+AW 



AA k v l 

8 (Av k ) 
dx» 



(C19) 



for this type of contraction. 

For the next contraction we return to the first split we 
used on the right-hand side: 

AWVpUu = [{n v + v v ) d p A + A V ^ + A V M ?v] , 

(C20) 

where again A 1 is spacelike. The first term gives simply 
AW[{n v + v v )d ll A] = -A i ^. (C21) 



From Eq. (jC8|) . the second term vanishes. In the third 
term we use Eqs. (|C9[) and (|A3|) . whence 



AW (AV M »„) = -AAWV„n v (C22) 
= AA^Ku. (C23) 



(Cll) Altogether, we have 



AWV u u v = AA l v k K kl - A 



dA 

dx l 



(C24) 



for this type of contraction. 

For the final contraction, one with a spacelike tensor 
B" 1 *, we return to the alternative 

B^V^Uu = B Ufl [n„ d p A + A V„n„ + V M (Av v )] (C25) 

on the right-hand side. The first term vanishes because 
B^^ is spacelike. The second term immediately yields, 
via Eq. (fA"3"]) . 



B"" (A V M n„) - AB kl V ink = -AB kl K kl . (C26) 
In the third term it once again will turn out best to raise 
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the index on v v . We find 

B"" [V M (A^)] - B^Vp (A«") (C27) 
= B„" [9 M (A^)+T^ (AO](C28) 
= 5**$ (A« fc ) 

+AB fcl « m .9^r^ (C29) 
= (A v fc ) 

+AB ki v m (Oakm 

+d m iki - dk ) (C30) 
= B k l d t (Av k ) 

+ AB ki v m d m7ki , (C31) 
where in the last step, the first and third terms — 



antisymmetric in i and k — vanish upon contraction with 
the symmetric B k% . Altogether, we have 



dx r < 



AB ki K, 



hi 



Bt 



; d(Av k ) 



(C32) 

for this type of contraction. 

In summary, Eqs. (fUTTj) . (JCTSJ, ([CM]) , and (jU32)) ex- 
hibit the four types of contractions appearing when either 
Eq. (|C3]) or (|U4|) is contracted with V^u,,. 
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